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SELF-MAGNETIC  FIELD  EFFECTS  ON  ELECTRON  EMISSION 
AS  THE  CRITICAL  IS  APPROACHED 


Abstract 

The  self-magnetic  field  associated  with  the  current  in  a  planar  diode  is  shown  to 
reduce  electron  emission  below  the  Child-Langmuir  current  density.  As  the  magnetic 
field  increases,  the  diode  current  is  limited  to  the  critical  current.  Here,  a  1-D  analysis  is 
carried  out  to  calculate  the  suppressed  current  density  in  the  presence  of  a  transverse 
magnetic  field.  The  problem  is  shown  to  be  similar  to  that  of  the  limiting  current  (i.e., 
Hull  current)  calculated  in  a  crossed  field  gap,  in  which  a  constant  transverse  magnetic 
field  is  applied  across  the  gap  to  insulate  the  electron  flow.  In  the  case  considered 
here,  the  magnetic  field  is  produced  by  the  diode  current  itself  and  this  self-magnetic 
field  decreases  with  distance  along  the  gap.  It  is  shown  that  the  emitted  current  density 
is  only  modestly  reduced  from  the  Child-Langmuir  current  density.  The  1-D  analysis 
remains  valid  until  critical  current  is  approached,  at  which  point  orbit  crossing  occurs 
and  a  2-D  kinetic  analysis  is  required.  The  minimum  diode  length  required  to  reach 
critical  current  is  also  derived. 


I.  Introduction 

High  power  vacuum  diodes  are  used  to  produce  intense  electron  beams  for  many 
applications.  When  self-magnetic  field  effects  are  negligible  and  the  cathode  is  a 
space-charge-limited  (SCL)  emitter,  Child-Langmuir  flow  is  obtained  [1,2].  Here, 
cylindrical  diodes  are  being  considered  with  an  anode-cathode  gap  that  is  much  less 
than  the  diode  radius.  In  this  case,  the  problem  can  be  treated  as  planar  and,  when  the 
magnetic  field  is  negligible,  electrons  cross  the  diode  with  straight-line  trajectories  as 
illustrated  in  Fig.  la.  However,  when  the  self-magnetic  field  associated  with  the  diode 
current  becomes  large  enough  to  significantly  deflect  the  electron  path  as  it  crosses  the 
diode  gap,  the  additional  electron  space  charge  density  in  the  gap  above  the  emission 
site  reduces  the  emission  below  the  Child-Langmuir  current  density.  This  situation  is 
illustrated  in  Fig.  1b.  At  high  magnetic  field,  the  diode  current  is  limited  to  the  critical 
current  ICnt,  which  is  obtained  when  the  electrons  reach  the  anode  at  grazing  incidence 
(or  equivalently  the  “electron  gyroradius”  equals  the  gap  size)  [3,4],  Critical  current  flow 
is  illustrated  in  Fig.  1c. 

Here,  a  1-D  analysis  is  carried  out  to  calculate  the  suppressed  current  density  in  the 
presence  of  a  transverse  magnetic  field.  The  problem  is  shown  to  be  similar  to  that  of 
the  limiting  current  (i.e.,  Hull  current)  calculated  in  a  crossed  field  gap  [5],  in  which  a 
constant  transverse  magnetic  field  is  applied  across  the  gap  to  insulate  the  electron 
flow.  In  treating  electron  flow  in  crossed  field  gaps,  the  Llewellyn  approach  is  typically 

Manuscript  approved  August  1 , 200 1 . 


1 


used  [5],  whereas  here,  Poisson’s 
equation  is  solved  by  direct  integration. 
Also,  in  the  case  considered  here,  the 
magnetic  field  is  produced  by  the  diode 
current  itself  and  this  self-magnetic  field 
increases  with  distance  along  the  gap 
causing  the  curvature  of  the  electron 
trajectories  to  increase  with  distance  as 
well  (see  Fig.  1c).  It  is  shown  that  the 
emitted  current  density  is  only  modestly 
reduced  from  the  Child-Langmuir  current 
density.  The  1-D  analysis  remains  valid 
until  critical  current  is  approached,  at 
which  point  orbit  crossing  occurs  and  a  2- 
D  kinetic  analysis  is  required. 

This  problem  is  of  interest  for 
determining  the  minimum  length  of  a  diode 
designed  to  run  at  lent-  This  length  will  be 
defined  as  the  critical  length  Lcrjt.  This 
situation  arises  when  an  electron  beam 
pinch  is  desired  to  concentrate  the 
electron-beam  energy,  while  minimizing 
losses  to  ion  current.  If  anode  plasma 
forms,  ion  current  is  generated  and  its 
magnitude  is  proportional  to  the  diode 
length  in  cylindrical  diodes  (or  diode  radius 
in  pinched-beam  diodes)[6-10].  Thus, 
minimizing  the  diode  length  minimizes  the 
current  lost  to  ions.  Although  this  analysis 
does  not  include  ion  current,  a  reasonable 
estimate  of  the  desired  minimum  diode 
length  with  ion  current  is  easily  obtained. 
For  a  nonrelativistic  planar  diode,  it  is  well 
known  that  the  SCL  electron  current 
density  is  enhanced  by  a  constant  factor  of 
1.86  when  ion  current  is  included  in  the 
analysis[1 1],  The  scaling  of  the  current 
density  with  diode  voltage  and  gap  are  the 


Figure  1.  Schematic  of  electron  flow 
accelerated  by  the  electric  field  E  across 
a  diode  with  a  transverse  self-magnetic 
field  B  that  is  generated  by  the  diode 
current.  Here,  V  is  the  diode  voltage,  Id0 
is  the  total  diode  current,  D  is  the  diode 
anode-cathode  gap  spacing,  and  L  is  the 
length  of  the  diode  emission  area.  The 
thin  arrows  indicate  the  electron 
trajectories  (a)  for  Id0  «  lent,  (b)  for  Id0  < 
lent,  and  (c)  for  Id0  =  Icrit-  A(z)  is  the  axial 
distance  an  electron  is  deflected  by  the 
local  transverse  magnetic  field  while 
crossing  the  gap. 

same  in  bipolar  flow  with  ions  as  it  is  when 


only  electron  current  flows.  Thus,  it  is  assumed  here  that  the  minimum  diode  length  for 


bipolar  flow  LB^t  is  to  first  order  simply  a  factor  of  1.86  smaller  than  the  minimum  length 


without  ions  (i.e.,  L^,  =  Lcrjt/1.86). 


In  Sec.  II  of  this  paper,  a  1-D  analysis  is  carried  out  to  calculate  the  electron  current 
density  as  a  function  of  the  transverse  magnetic  field.  The  results  of  this  analysis  are 
then  applied  in  Sec.  Ill  to  find  the  current  as  a  function  of  distance  along  a  diode.  By 
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equating  this  current  to  Icrit,  the  critical  length  is  obtained.  Finally,  the  results  of  this 
work  are  summarized  in  Sec.  IV. 


II.  Electron  current  density  as  a  function  of  the  transverse  magnetic  field 

The  planar  diode  geometry  considered  here  is  illustrated  in  Fig.1.  A  voltage  V  is 
applied  across  the  electrodes  that  are  separated  by  a  distance  D.  For  this  part  of  the 
analysis,  a  uniform  magnetic  field  is  applied  transverse  to  the  electron  flow  in  the  y 
direction  and  current  is  emitted  over  an  axial  length  L  as  shown  in  Fig.  la.  For  this 

situation,  the  electric  potential  O(x)  and  the  magnetic  vector  potential  Az(x)  can  be 
written  as 

d°(x)  -  E  M 

dx  Ex(  }  '  0) 

Az(x)  =  -B0x  ,  ^2) 

where  Ex  is  the  electric  field  in  the  gap,  the  constant  magnetic  field  B0  >  0,  0(0)  =  0,  and 
O(D)  =  V.  For  SCL  flow,  an  additional  boundary  condition  is  Ex(0)  =  dO/dx(0)  =  0.  In 
equilibrium,  the  continuity  equation  requires  that  the  current  density  Jx(x)  be  given  by 

Jx(x)  =  -Jo  .  (3) 

where  again  the  constant  Jo  >  0.  Because  y  and  z  are  ignorable  coordinates,  the 
canonical  momentum  associated  with  each  of  these  coordinates  is  conserved.  Thus, 
the  velocity  components  vy(x)  and  vz(x)  of  the  electrons  can  be  written  as 
vy(x)  =  0  , 

vz(x)  =  -©cex  ,  (5) 

where  coce  =  eB0/mec,  e  is  the  electron  charge,  me  is  the  electron  mass,  c  is  the  speed  of 
light,  and  it  is  assumed  that  the  electrons  are  born  on  the  cathode  surface  at  x  =  0  with 
zero  energy.  Using  this  same  assumption,  an  expression  for  vx(x)  is  obtained  from 
conservation  of  energy, 

,(X)  =  2e4>M 


-<x2 


Here,  Eqs.  (4)  and  (5)  have  been  used  to  substitute  for  vy  and  vz. 

Equation  (6)  immediately  provides  an  expression  for  the  critical  current  Icrit-  Because 
vx  >  0  is  required  for  a  physically  meaningful  result  from  Eq.  (6),  the  current  I  must  be 
less  than  or  equal  to  Icrit,  where 


^ciit 


R  -  21 

Bo  =  f^ 


Note  that  B0  is  defined  as  in  Eq.  (8)  so  that  these  ID  results  can  be  related  to  problems 
in  cylindrical  geometry  with  R  »  D,  where  R  is  the  cathode  radius.  From  Eqs.  (7)  and 
(8),  the  critical  magnetic  field  is  found  to  be  Bocrit  =  2Icrjt/Rc.  The  critical  current  is  also 
referred  to  as  the  Hull  current  for  systems  where  Bo  is  applied  [5]. 
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An  equation  for  O(x)  can  now  be  obtained, 
electron  density  as 


ne(x)  = 


evx(x) 

Poisson’s  equation  becomes 
d20  4nJn 


dx5 


2eO 

itl 


2  2 
~®ceX 


Using  Eq.  (3)  and  Jx  =  -enevx  to  write  the 


(9) 

(10) 


The  solution  to  this  equation  that  satisfies  the  boundary  conditions  (0(0)  =0,  O(D)  =  V, 
and  dO/dx(0)  =0)  provides  a  complete  solution  to  the  problem  by  substituting  O(x)  back 
into  Eqs.  (1),  (6),  and  (9)  for  Ex(x),  vx(x),  and  ne(x)  respectively.  In  addition,  vy(x)  and 
vz(x)  are  given  by  Eqs.  (4)  and  (5).  Because  Eq.  (10)  is  a  second-order  differential 
equation,  it  requires  only  two  boundary  conditions.  The  third  boundary  condition 
determines  the  eigenvalue  for  J0. 


In  order  to  solve  Eq.  (10),  it  is  convenient  to  rewrite  the  equation  as 
d2S  _  87teJ0  0  2 

where 


(11) 


s  =  vx  ■  (12) 

and  vx2  is  given  in  terms  of  O(x)  and  x  in  Eq.  (6).  The  first  integral  can  be  obtained  by 

multiplying  both  sides  of  the  equation  by  dS/dx  and  integrating.  This  yields 

\  1/2 


dS 

dx 


=  +! 


^^Vs-teLs 

m 


(13) 


where  S(0)  =  0  and  dS/dx(0)  =  0  have  been  used.  These  new  boundary  conditions  can 
be  derived  from  the  boundary  conditions  on  O  at  x  =  0,  and  Eqs.  (12)  and  (6).  Because 
S  is  associated  with  vx2  [see  Eq.  (12)],  the  upper  sign  applies  in  the  usual  case  when  vx 
is  increasing,  while  the  lower  sign  applies  when  vx  is  decreasing.  Note  that  vx  will  only 
decrease  if  the  electron  orbit  has  been  significantly  deflected  by  the  magnetic  field  so 
that  energy  is  being  transferred  into  axial  motion.  This  occurs  when  I  approaches  Icrit- 


The  second  integral  can  be  obtained  by  directly  integrating  Eq.  (13)  and  using  the 
boundary  condition  S(0)  =  0.  Using  Eqs.  (6)  and  (12)  to  write  the  solution  in  terms  of 
O(x),  yields  for  dS/dx  >  0  [i.e.,  the  upper  sign  in  Eq.  (13)] 


x  =  — 


87teJ0 1 
me®ce  l 

2eO(x)  2 

m  ®ceX 

me 

47ieJ0 

Sin"1 

me®ce  ( 

me«ce  * 

47teJ0  v 

.1/2 


O), 


^2eO(x)  , 

— T 

V  me  J 


1/2 


2eO(x) 


,1/2 


rri 


2  2 
-“ceX 


-1 


71 

+  — 


(14) 
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This  solution  only  applies  for  x  <  x(1)  where  x(1)  is  defined  below  [see  Eq.  (24)]  as  the 
position  where  vx  stops  increasing  (i.e.,  where  dS/dx  =  0).  The  value  of  x(1)  depends  on 
the  magnetic  field  strength  through  ©<*.  When  x(1)  >  D,  Eq.  (14)  applies  for  all  0  <  x  <  D. 
When  x(1)  <  D,  Eq.  (14)  applies  for  0  <  x  <  x(1)  and  Eq.  (13)  must  be  solved  again  using 
the  lower  sign  this  time  for  x(1)  <  x  <  D  and  that  solution  matched  to  the  solution  in  Eq. 
(14)  at  x  =  x  .  This  solution  is  found  below  [see  Eq.  (26)].  Assuming  for  now  that  x^  > 
D,  the  boundary  condition  <h(D)  =  V  can  be  applied  to  Eq.  (14)  to  provide  a 
transcendental  equation  for  the  eigenvalue  J0  for  a  given  value  of  a^.  Once  Jo(©ce)  is 
determined,  Eq.  (14)  can  be  solved  numerically  for  O(x). 


It  is  instructive  to  examine  this  solution  for  small  values  of  B0.  To  accomplish  this,  let 

m  yv.2  f  Ort,f././\ 


47ieJn 


2eO(x) 


m« 


2  2 


«1 


Using  the  expansion 


,3/2 


7s 


.5/2 


Sin-’<e  - 1)  =  -  +  (2e)  "•  +  — T t 
2  6V2  144V2 


+  .... 


for  e  « 1 


Eq.  (14)  can  be  expanded  to  give 
/  \1/2  / 

m.  I  1 2eO(x) 

l  m. 


x  = 


18rceJ 


-®ceX2 


0  J 


\  3/4 

2  f 

1  ,  me®ce 

) 

247ieJ0  v 

2eO(x) 


V'2 


m. 


(15) 


(16) 


(17) 


First,  note  that,  for  B0  =  0  (i.e.,  ace  =  0),  the  Child-Langmuir  equation  is  recovered.  To 
obtain  this  result,  the  final  boundary  condition,  O(D)  =  V,  is  applied  yielding 


JS”=Jc8-SL 

9tc 


(  \1/2  y3/2 


2e 

vme; 


D2 


(18) 


Here,  is  defined  as  the  value  for  Jo  at  ©ce  =  0.  Equation  (17)  can  also  be 

manipulated  to  find  an  iterative  solution  for  J0  with  x  =  D  and  <D(D)  =  V.  Retaining  only 
the  first  two  terms  gives 

f  2  nv2 

rmaiD2 ] 


Jn  = 


971 


2e 


V-' 


2e 


D" 


l4  3me©^D2/4e 


v  _  me©^D 


2\ 


2e 


(19) 


This  expansion  is  the  Child-Langmuir  current  corrected  for  the  presence  of  the  small 
magnetic  field,  which,  as  expected,  reduces  J0  below  JCL.  Although  the  leading  terms  in 
the  expansion  shown  in  Eq.  (19)  vanish  as  I  ->•  Icrit,  the  higher  order  terms  in  the 
expansion  that  are  not  shown  diverge  as  I  -»  Icrit  because 

2  ‘  (20) 


.  //..crit\2  _  2eV 

“ce  “>  (“ce  )  = - 7^2 

m„D 


Thus,  Eq.  (19)  is  only  valid  for  ©ce  <  © 


crit 

ce 
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As  the  magnetic  field  strength  increases,  a  point  is  reached  where  the  electron  orbit 
has  been  deflected  axially  so  much  that  vx  begins  to  decrease  before  the  electron 
reaches  the  anode.  In  fact,  for  ©ce  =  ©£f ,  the  electron  slows  down  in  the  x  direction  to 

the  point  that  vx  =  0  when  the  electron  reaches  the  anode.  When  vx  stops  increasing, 
dS/dx  =  0,  and  when  vx  begins  to  decrease,  the  lower  sign  in  Eq.  (13)  applies.  To  find 
the  magnetic  field  strength  where  the  electron  reaches  the  anode  just  as  the  electron 
stops  accelerating  the  x  direction,  Eq.  (13)  is  solved  at  x  =  D  (and  O(D)  =  V)  with  dS/dx 
=  0.  This  yields 


8weJW 


m. 


=«’F 


(21) 

where  g>™  denotes  the  magnetic  field  strength  where  dS/dx  =  0  at  x  =  D  and  J™  is  the 
current  density  for  this  special  magnetic  field  strength.  This  result  can  be  substituted 
into  Eq.  (14)  to  obtain  an  expression  for  J™  with  ©«  =  ©™ .  Substituting  this  result  back 
into  Eq.  (21)  and  solving  for  ©™  yields  for  x(1)  =  D 


and 


2  eV 
vmeD2, 


,1/2 


0.844©^ 


(22) 


meD 


(<)3=- 


9tc2J 


CL 


4K2e'”“;  '2(n2  +  4)”!  ■3°'86Jcl  '  <23> 

Thus,  for  0  <  ©ce  <  ©^ ,  Eq.  (14)  applies  and  J0  decreases  from  JCl  at  ©<*  =  0  to  0.86JCl 
at  ©ce  =  (0™  • 


For  ©™  <  ©ce  ^  ©«  ,  the  problem  is  more  complicated.  Equation  (14)  applies  for  0  < 

x  <  x(1)  where  x(1)  is  defined  to  be  the  position  where  dS/dx  =  0  (i.e.,  where  vx  begins  to 
decrease).  However,  for  x(1)  <  x  <  D,  Eq.  (13)  must  be  solved  again  using  the  lower  sign 
and  then  this  new  solution  must  be  matched  to  Eq.  (14)  at  x  =  x^\  An  expression  for 
x(1)  can  be  found  by  equating  the  right  hand  side  (RHS)  of  Eq.  (13)  to  zero  and 
substituting  the  resulting  value  of  S(x(1))  into  Eq.  (14)  evaluated  at  x  =  x(1).  This 
procedure  yields 


x(1)  = 


4n2eJ0 


(24) 


and 


S(x(1>)  = 


8neJ0 
me«ce  > 


(25) 


where  J0  is  a  function  of  ©ce  and  is  yet  to  be  determined.  The  matching  condition  at  x(1) 
can  then  be  used  to  give 
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x  =  + 


8rceJn 


me< 


4rceJ0 

me®L 


2eO(x) 


J/2 


-®Lx2 


© 


2eO(x)  2  2 


m. 


1/2 


Sin 


-1 


47ieJr 


2eO(x) 
m„ 


-<x2 


1/2 


-1 


371 


(26) 


V  ) 

for  x  <  x  ^  D  and  ©2  <  ®ce  ^  ©«  •  Applying  the  boundary  condition  O(D)  =  V  at  x  =  D 
to  Eq.  (26)  provides  a  transcendental  equation  for  the  eigenvalue  J0  for  a  given  value  of 
coce-  Note  that,  for  the  special  case  of  ©«*  =  ©  J ,  Jo  is  defined  as  Jf  which,  from  Eqs. 
(18),  (20),  and  (26),  is  found  to  be 

9J“  ->0.716Jcl  .  (27) 

Once  Jo  is  determined  as  a  function  of  ©<*,  Eqs.  (14)  and  (26)  can  be  solved  numerically 
for  O(x)  as  a  function  of  ©ce. 


icrit  _ 

0  4k 


In  order  to  examine  the  behavior  of  the  solution  near  ©ce  =  ©£f ,  Eq.  (26)  can  be 
expand  about  x  =  D  using  the  same  small  parameter  e  as  defined  in  Eq.  (15).  Here,  e  is 
small  because  2eO/me  ~  ©^  x2  for  ©ce  ~  ©£f  and  x  ~  D.  As  done  with  J0  for  ©ce  ~  0  in 
EQ-  (19),  this  expansion  can  also  be  manipulated  to  find  an  iterative  solution  for  J0  with  x 
=  D  and  O(D)  =  V  for  ©ce  ~  ©”f .  Retaining  only  the  first  two  terms  gives 


Jn  = 


9J 


CL 


4n 


1  +  -f— 

3  2 


2  V'V 


1- 


© 


ce 


CO 


crit 
ce  J 


-3 


1- 


®ce 

-f .... 

^Crit 

®ce 

J 

where  Eq.  (27)  is  recovered  for  ©ce  =  ©^ .  Note  that,  as  ©«,  -»  ©* 
with  respect  to  ©ce  diverges  as 


aj0 

_  ^CL 

r 7i 2 

5©ce 

5  8 
3 

K 

tJ- 

1 

U 

.  1/4 


f  V'4 

\  ®ce 
crit 

’ce  J 


+  ....  — >  —oo 


OUJce  ^®ce  V  *  J  \  ®o e  J 

Thus,  the  electron  emission  abruptly  shuts  off  at  ©ce  =  © 


-  „crit 
ce 


(28) 

the  derivative  of  J0 
(29) 


A  numerical  solution  of  the  transcendental  equation  for  J0  has  been  obtained  for  all 
values  of  ©ce  in  the  range  0  <  ©ce  <  ©*‘ .  For  0  <©<*<©£ ,  the  transcendental  equation 
given  in  Eq.  (14)  is  used  with  O  =  V  and  x  =  D,  while,  for  ©£><  ©ce  <  ©^ ,  Eq.  (26)  is 

used  with  the  same  boundary  conditions.  A  plot  of  Jo  as  a  function  of  ©ce  is  shown  in 
Fig.  2  as  the  solid  curve.  As  expected  from  the  above  results,  the  current  density 
decreases  slowly  with  magnetic  field  until  ©^‘  is  reached  where  emission  abruptly  shuts 

off.  The  approximation  for  J0  given  in  Eq.  (19)  is  also  plotted  in  Fig.  2  as  the  dashed 
curve.  This  approximation  is  only  reasonably  accurate  for  <  0.7.  Retaining  additional 
terms  in  the  expansion  would  extend  the  accuracy  of  the  approximation  to  larger 
©ce/oC  but  the  additional  terms  would  diverge  at  ©ce/©^  =  1-  Note  that  retaining  only 
the  first  term  in  the  expansion  given  in  Eq.  (19),  provides  an  approximation  that  is  only 
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accurate  for  very  small  ©ce/©^1  (see 
dotted  curve  in  Fig.  2).  A  plot  of 
x(1)/D  as  a  function  of  ©ce  is  also 
plotted  in  Fig.  2  where  x(1)  <  D. 

Using  the  solution  for  J0  shown 
in  Fig.  2,  Eqs.  (14)  and  (26)  can  be 
used  to  solve  for  d>(x)  for  a  given 
value  of  ©ce-  Equations  (6),  (5),  and 
(9)  can  then  be  used  to  obtain  vx(x), 
vz(x),  and  ne(x),  respectively  for  a 
given  value  of  ©«,.  Solutions  for 
<D(x),  ne(x),  vx(x),  and  vz(x)  are 
plotted  in  Figs.  3-6  respectively  for 
cases  with  ©<*/©«  =  0,  0.25,  0.5, 

0.75,  and  1.0.  In  these  plots,  vx(x), 
and  vz(x)  are  normalized  to  v0  = 
(2eV/me)1/2,  and  ne(x)  is  normalized 
to  no  =  ne(D)  for  ©ce/ ©^l  =  0.  From 

Eq.  (18),  n0  =  Jci_/ev0  =V/(9rceD2). 
The  profile  of  the  potential  4>(x) 
does  not  vary  much  as  ©ce/©« 


Figure  2.  Plots  of  the  normalized  electron 
current  density  Jo/Jcl  (solid)  and  the  normalized 
distance  x(1)/D  for  x(1)  <  D  (dot-dash)  as 
functions  of  the  normalized  transverse  magnetic 
field  strength  ©ce/®cecrit  =  B0/B0crit.  Also  shown 
are  plots  of  J(/0)  (dot),  the  zero  order 
approximation  to  J0  from  the  first  term  in  Eq. 
(19),  and  J0(1)  (dash),  the  first  order 
approximation  to  J0  from  the  first  two  terms  in 
Eq.  (19). 


Figure  3.  Plots  of  the  normalized 
potential  0(x)/V  as  a  function  of  the 
normalized  distance  across  the  gap  x/D 
for  various  values  of  ©ce/fflcecrit-  Plots  for 
®ce/®cecnt  =  0.0  (solid),  0.25  (dash),  0.5 
(dash-dot),  and  0.75  (dash-dot-dot)  nearly 
lay  on  top  of  each  other,  while  the  plot  for 
®ce/®cecnt  =1.0  (dot)  lies  slightly  below  the 
other  curves. 


0.0  0.2  0.4  0.6  0.8  1.0 


x/D 

Figure  4.  Plots  of  the  normalized  density 
ne/n0  as  a  function  of  the  normalized 
distance  across  the  gap  x/D  for  various 
values  of  ©ce/©ceCnt.  Plots  for  ©ce/©ceCrit  = 
0.0  (solid),  0.25  (dash),  0.5  (dash-dot), 
and  0.75  (dash-dot-dot)  nearly  lay  on  top 
of  each  other,  while  the  plot  for  wce/®cecrit  = 
1.0  (dot)  diverges  as  it  approaches  the 
anode  at  x  =  D. 
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Figure  5.  Plots  of  the  normalized  velocity 
vx/v o  as  a  function  of  the  normalized 
distance  across  the  gap  x/D  for  various 
values  of  coce/©cecnt.  Plots  are  shown  for 
©ce/ooce0"1  =  0.0  (solid),  0.25  (dash),  0.5 
(dash-dot),  0.75  (dash-dot-dot),  and  1.0 
(dot). 


Figure  6.  Plots  of  the  normalized  velocity 
vz/v0  as  a  function  of  the  normalized 
distance  across  the  gap  x/D  for  various 
values  of  coce/cocecrit.  Plots  are  shown  for 
©ce/©cecnt  =  0.0  (solid),  0.25  (dash),  0.5 
(dash-dot),  0.75  (dash-dot-dot),  and  1.0 
(dot). 


increases  with  four  of  the  five  curves  lying  nearly  on  top  of  each  other  and  the  coce/©S‘  = 
1  curve  lying  only  slightly  below  the  other  curves.  The  density  profile  ne(x)  is  also  little 
changed  over  most  of  the  range  for  ©ce/©^‘  but  changes  dramatically  as  ©ce/©^‘ 
approaches  unity  where  the  density  diverges  at  the  anode.  The  profiles  of  the  electron 
velocity  components  vx(x)  and  vz(x)  both  depend  strongly  on  gWcd^.  For  ©<*/©£?  =  1, 
note  that  vx(x)  goes  to  zero  as  the  electrons  approach  the  anode  at  grazing  incidence. 


III.  Minimum  diode  length  for  obtaining  critical  current 

If  laminar  flow  is  assumed  in  the  diode  gap,  then  the  above  analysis  can  be  used  to 
calculate  the  diode  current  Id(z)  as  a  function  of  distance  along  the  cathode.  The  total 
current  in  the  diode  is  Id0  as  shown  in  Fig.  1 .  Id(z)  increases  from  zero  at  z  =  0  to  Id0  at  z 
=  L.  Following  the  notation  used  in  Eq.  (8),  the  local  self-magnetic  field  associated  with 
this  diode  current  is  given  by 

B°<z,=i*r  ’  (30) 

and  determines  the  local  current  density  from  J0  as  a  function  of  B0  given  in  Fig.  2. 
Here,  this  local  current  density  is  denoted  by  J0(z)  and  dld(z)/dz  =  2nRJo(z).  This  diode 
configuration  is  illustrated  in  Fig.  la  for  Id0  «  Icrit,  in  Fig.  1b  for  Id0  <  Icrit,  and  in  Fig.  1c 
for  Id0  =  Icrit-  This  analysis,  which  assumes  that  z  is  an  ignorable  coordinate,  is  still 
approximately  correct  as  long  as  variations  in  the  z  direction  are  small  compared  to 
variations  in  the  x  direction  (i.e.,  d/dz  «  d/dx).  This  is  supported  by  the  fact  that  J0(z) 
only  varies  slowly  with  Bo(z)  as  shown  in  Fig.  2. 


9 


The  minimum  diode  length  needed  to  obtain  critical  current  then  is  determined  by 
Id(Lcrit)  =  2kR^c,#  J0(z)dz  =  ICTit  ,  (31) 

where  ICrit  is  defined  in  Eq.  (7).  Because  Jo  decreases  slowly  with  ©ce/co^l*  (i.e.,  with  B0), 

a  zero  order  estimate  (and  lower  bound)  for  Unt  can  be  obtained  by  assuming  that  J0(z) 
is  uniform  in  z  and  equal  to  JCl,  yielding 
,0  _  9D 

crit  2p2  '  (32) 

Here,  p  =  Vo/c.  Note  that  L°crit  »  D  for  nonrelativistic  beams,  and  that  L°crit  ->•  9D/2  as 
P^I  where  this  nonrelativistic  analysis  breaks  down.  Using  L0crit  to  scale  z,  Eq.  (31) 
can  be  integrated  to  solve  for  Lcrit.  The  results  are  plotted  in  Fig.  7  and  show  that 
Lcm  =1.076^,  .  (33) 

This  provides  both  scaling  for  and  a  reasonable  estimate  of  the  minimum  diode  length 
required  for  critical  current. 


To  examine  the  validity  of  the  laminar  flow  assumption,  consider  the  distance  A(z) 
that  an  electron  is  deflected  axially  while  traversing  the  diode  gap  as  illustrated  in  Fig. 
1b.  At  z  =  0,  the  self-magnetic  field  vanishes  so  A(0)  =  0.  For  z  >  0,  A(z)  monotonically 
increases  and  can  be  written  as 


A(z)  =  £  vz(x)dt  =  £  ^dx 


where  T  is  the  transit  time  of  the  electron  for  crossing  the  diode  gap  D,  and  vz  and  vx  are 
defined  in  Eqs.  (5)  and  (6)  respectively.  Noting  that  the  denominator  of  the  integrand  is 
just  vx  =  S1/2  [see  Eq.  (12)]  and  that  Eqs.  (14)  and  (26)  can  be  used  to  write  vz(x)  =  ©ceX 
in  terms  of  S,  Eq.  (34)  can  be  rewritten  as  an  integral  in  S  where  from  Eq.  (13) 

dx  =  +^^eJo  Vs-4©^s  dS  .  (35) 

When  0  <  ©ce  <  ©^ ,  Eq.  (14)  is  used  to  substitute  into  Eq.  (34)  for  x  to  write  the  integral 

in  terms  of  S.  When  ©^  <  ©ce  ^  ,  the  integral  must  be  broken  into  two  parts;  for  0 

<  x  <  x(1),  Eq.  (14)  is  used  to  substitute  for  x  in  terms  of  S,  and  for  x(1)  <  x  <  D,  Eq.  (26)  is 
used.  Performing  the  integration  yields 


A(z)  = 


1  ( 2eV 


-©2eD; 


2rc2eJ„ 


me©ce 


me©ce 


me©ce 


—  +  2]  Sin’1 
2  47teJA 


for  0  <  ©ce  <  ©(^ ,  and 
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Figure  7.  Plots  of  the  normalized  current 
density  J0(z)/JCL  (solid)  and  the 
normalized  diode  current  Id(z)/Icrit  (dash) 
as  functions  of  the  normalized  distance 
along  the  diode  z/L°crit. 


Figure  8.  Plots  of  the  normalized  axial 
deflection  of  an  electron  A(z)/D  (solid)  as  it 
crosses  the  gap  and  the  normalized 
derivative  of  the  axial  deflection 
d(A(z)/D)/d(z/L0crit)  (dash)  as  functions  of  the 
normalized  distance  along  the  diode  z/L0^. 


-if  2eV  2r,2 
A(z)  = - - ©*  D2 


^eJo  o:„-l  me©L  ( 2eV 


47teJ0  ^  me 


©LD2  -l 


^4-f  — + 2  j  Sin"1  [  2 eV  V" 

me®L  2  4rceJ0  ^  me  06  J 


for  ©^  <  ©ce  <  .  In  Eqs.  (36)  and  (37),  the  z  dependence  appears  on  the  RHS 

through  J0(z)  and  ©ce(z)  =  eB0(z)/mec,  where  the  explicit  dependence  on  z  has  been 
suppressed  for  the  sake  of  brevity.  Note  that  Eqs.  (36)  and  (37)  can  be  written  in 
dimensionless  form  by  using  the  normalized  parameters  A(z)/D,  J0/JCl,  and  ©ce(z)/ffi^ . 


The  orbit  of  an  electron  emitted  from  the  cathode  at  z  +  6  crosses  the  orbit  of  an 

fKeCoP,oeTi!!ed  at  z  lf  z "  A(z)  >  (z  +  8)  -  A(z  +  6).  In  the  limit  of  8  «  z,  the  last  term  on 
the  RHS  of  this  equation  can  be  Taylor  expanded.  Upon  substitution  of  this  expansion, 

the  leading  terms  in  the  orbit  crossing  condition  cancel.  Thus,  in  the  limit  of  8  ->  0,  the 
criteria  for  orbit  crossing  simply  becomes 
dA 

dl  •  (38) 

In  order  to  evaluate  Eq.  (38),  expressions  for  dJ0/dz  and  d©ce/dz  are  needed.  First  note 
that  dJ0/dz  -  (dJ0/d©ce)d©ce/dz  where  the  solution  for  Jo(©ce)  is  shown  in  Fig.  2  and  its 
differentiation  with  respect  to  ©ce  is  straight  forward.  Then  note  that  d©ce/dz  = 
(2e/meRc  )dld/dz  =  4neJ0(z)/mec2  where  Id(z)  and  J0(z)  are  shown  in  Fig.  7.  Plots  of  A(z) 


11 


and  dA(z)/dz  are  shown  in  Fig.  8  where  A(z)  has  been  normalized  to  D 
normalized  to  L°Crit-  With  this  normalization  Eq.  (38)  becomes 
d(A(z)/D)  L°crit_  9 
dCz/O  D  2P2 


and  z  has  been 


(39) 


Thus,  as  the  voltage  increases  (where  p2  ~  V),  the  location  where  orbit  crossing  first 
occurs  moves  to  smaller  z.  The  case  with  p  =  0.5  (i.e.,  V  =  64  kV)  yields  L°crit  /D  =  18 
and  is  shown  as  an  example  in  Fig.  8.  For  this  case,  orbit  crossing  first  occurs  at  z  ~ 
1 .055L°crit.  Because  orbit  crossing  only  occurs  very  close  to  Unt  =  1  076  1%,  the  ana 
lysis  is  valid  over  nearly  the  full  diode  length  and  Eq.  (33)  provides  an  accurate  estimate 
for  Lent- 


IV.  Conclusions 

The  self-magnetic  field  associated  with  the  diode  current  in  a  planar  diode  is  shown 
to  reduce  electron  emission  below  the  Child-Langmuir  current  density.  As  the  magnetic 
field  increases,  the  diode  current  is  limited  to  the  critical  current  Icrit-  This  problem  is  of 
interest  for  determining  the  minimum  length  Lcrit  of  a  diode  designed  to  run  at  Icrit.  This 
situation  arises  when  an  electron  beam  pinch  is  desired  to  concentrate  the  electron- 
beam  energy,  while  minimizing  losses  to  ion  current.  If  anode  plasma  forms,  ion  current 
is  generated  and  its  magnitude  is  proportional  to  the  diode  length.  Thus,  minimizing  the 
diode  length  minimizes  the  current  lost  to  ions.  Here,  a  1-D  analysis  is  carried  out  to 
calculate  the  suppressed  electron  current  density  in  the  presence  of  a  transverse 
magnetic  field  without  ions.  The  problem  is  shown  to  be  similar  to  that  of  the  limiting 
current  (i.e.,  Hull  current)  calculated  in  a  crossed  field  gap,  in  which  a  constant 
transverse  magnetic  field  is  applied  across  the  gap  to  insulate  the  electron  flow.  In 
treating  electron  flow  in  crossed  field  gaps,  the  Llewellyn  approach  is  typically  used, 
whereas  here,  Poisson’s  equation  is  solved  by  direct  integration.  Also,  in  the  case 
considered  here,  the  magnetic  field  is  produced  by  the  diode  current  itself  and  this  self- 
magnetic  field  varies  with  distance  along  the  gap.  Special  values  of  some  of  the 
important  parameters  are  listed  in  Table  1. 


Table  1 .  Special  values  for  various  important  parameters.  Note  that 


JCL  =(2e/me)1/2V3/2 /97iD2,  v0  =  (2eV/me)1/2,  co^1  =v0/D,  and  n0=JCL/ev0. 
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It  is  shown  that  the  emitted  current  density  is  only  modestly  reduced  from  the  Child- 
Langmuir  current  density  JCl  with  the  current  density  only  dropping  to  0.71 6 JCl  at  1^. 
Application  of  the  1-D  analysis  to  the  pinched  beam  diode  remains  valid  until  critical 
current  is  approached,  at  which  time  orbit  crossing  occurs  and  a  2-D  kinetic  analysis  is 
required.  However,  it  is  shown  that  orbit  crossing  only  occurs  very  close  to  the  end  of 
the  pinched  beam  diode.  Thus,  a  reasonable  estimate  for  the  critical-current  diode 
length  is  Unt  =  1.0761.%  =  4.84D/p2.  Although  this  analysis  does  not  include  ion 
current,  an  argument  is  made  for  extending  these  results  for  the  case  with  ions.  For  a 
nonrelativistic  planar  diode,  it  is  well  known  that  the  SCL  electron  current  density  is 
enhanced  by  a  constant  factor  of  1.86  when  ion  current  is  included  in  the  analysis.  The 
scaling  of  the  current  density  with  diode  voltage  and  gap  are  the  same  in  bipolar  flow 
with  ions  as  it  is  when  only  electron  current  flows.  Thus,  it  is  assumed  here  that  the 
minimum  diode  length  for  bipolar  flow  L^t  is  to  first  order  simply  a  factor  of  1.86  smaller 

than  the  minimum  length  without  ions  (i.e.,  LB^  =  Lcrit/1.86  =  2.6D/p2).  This  assumes 

that  ion  current  is  much  less  than  the  electron  current.  Also  note  that  this  last  result 
only  applies  for  D  «  R,  because  it  has  been  shown  recently  that  the  factor  of  1.86 
increases  with  aspect  ratio  for  cylindrical  diodes  (with  the  cathode  radius  larger  than  the 
anode  radius)[12]. 

Extending  these  results  to  the  relativistic  regime  and  including  ions  in  the  analysis 
will  require  a  numerical  solution  of  Poisson’s  equation  and  is  the  subject  for  future  work. 
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